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Abstract 

The aim of this review is to present the Ust of by now a significant collection of quantum 
integrable models, ultralocal as well as nonultralocal, in a systematic way stressing on their un- 
derlying unifying algebraic structures. We restrict to quantum and statistical models belonging 
to trigonometric and rational classes with 2x2 Lax operators. The ultralocal models can be 
classified successfully through their associated quantum algebra and are governed by the Yang- 
Baxter equation, while the nonultralocal models, the theory of which is still in the developmental 
stage, allow systematization through the braided Yang-Baxter equation. Along with the known 
integrable models some possible directions for investigation in this field and generation of such 
new models are suggested. 

Key words: Integrable Quantum and Statistical Vertex Models; Quantum Algebras, Yang-Baxter 
Equation, Braided Extention, Algebraic Bethe Ansatz, Ultralocal and Nonultralocal Models 

1 Introduction 

By quantum integrable systems we will mean the systems with sufficient number of higher conserved 
quantities including the Hamiltonian of the model. Such a notion of integrability in the Liouville 
sense allows description through action-angle variables with the conserved quantities, which are now 
operators, playing the role of action variables. For integrable systems the conserved quantities, being 
functionally independent should form a commuting set of operators [c„ , Cm] =0, n, m = 1, 2, . . . , A^, 
such that their total number would match with the degree of freedom of the system. For example, 
an one dimensional lattice model of /-sites describing d-mode pseudo-particle should have number of 
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conserved quantities N = dl. Note that for spin-^ chains we have d = I, while spin-1 and electron 
models account for d = 2. In this review we will stick to single-mode: d = 1 systems for simplicity 
and consider mainly periodic lattice models with N < oo, where the algebraic structures can be seen 
in their exact form. At the lattice constant A — the field models will be generated from their exact 
lattice versions, whenever possible. Integrable field models with N ^ oo consequently needs to have 
infinite number of conservation laws. 

Integrable systems therefore are restrictive systems with a very rich symmetry. The beauty of such 
models is that they allow exact solutions for the eigenvalue problem simultaneously for all conserved 
operators including the Hamiltonian. Moreover such 1-dimensional quantum systems are related also 
to the corresponding 2-dimensional classical statistical models with a fluctuating variable. Therefore 
parallel to a quantum mechanical model one can in principle exactly solve also a related vertex-type 
model on a 2-dimensional lattice using almost the same techniques and similar results [1] . Celebrated 
examples of such interrelated integrable quantum and statistical systems are XYZ quantum spin-^ 
chain and the 8-vertex statistical model, XXZ spin chain and the 6-vertex model, spin-1 chain and 
the 19-vertex model etc. 

For describing an integrable system with such an involved structure, one naturally can no longer 
start from the Hamiltonian of the model as customary in physics, since now the Hamiltonian is merely 
one among many commuting conserved charges. It therefore needs to adopt certain abstractions 
which are formalized by the quantum inverse scattering method and the algebraic Bethe ansatz (see 
[2, 3]). Though we would use the same language, we take here a slightly different view point since 
we intend to describe integrable systems belonging to both ultralocal and nonultralocal classes. 
For effective description of integrable systems it is convenient to define a generating function called 
transfer matrix r(A), depending on some extra parameter A known as the spectral parameter, such 
that one can recover the infinite number of conserved quantities as the expansion coefficients of t(A) 
or any function of it like InT(A) = J^j Cj A-^. The crucial integrability condition may then be defined 
in a compact form as 

[r(A),r(^)]=0, (1.1) 

from which the commutativity of Cj's follows immediately by comparing the coefficients of different 
powers of A, n. 

However for solving the eigenvalue problem as well as for identifying the structure of the model we 
require a more general matrix formulation, from where the integrability condition may be derived. At 
the same time we need to transit from the global to the local description defined at each lattice point, 
where some individual properties of a model are well expressed. At this local level, as we see now, 
the difference between the ultralocal and the nonultralocal models become prominent. An integrable 
system allowing the needed abstraction may be represented by an unusual type of matrix called the 
Lax operator Laj (A) defined at each site j in a 1-dimensional discretized lattice. The index a defines 
the matrix or the auxiliary space, while j designates the quantum space. The matrix elements of the 
Lax operator, unlike in usual matrices are operators acting on some Hilbert space. The models with 
the Lax operators commuting at different lattice sites: 

LaM)Lbk{lA = Lbk{l^)Laj{X), a^b, j (1.2) 
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are known as the ultralocal models, while the integrable models for which the above ultralocality 
condition does not hold are classified as the nonultralocal models. Note that in expressions like 
(1.2) different auxiliary spaces mean different tensor products like Lij(X) = Lj(X) ® I and L2j(iJ,) = 
I®Lj{ij). The ultralocal property (1.2) generally reflects the involvement of canonical operators with 
commutation relations like [u{x)^p{y)\ = iS{x — y) or [■0(x), ■0^(y)] = S{x — y) in the Lax operator 
giving trivial commutator at points x ^ y. ha. nonultralocal models on the other hand the basic 
fields may be of noncanonical type, e.g. [jiix), ji{y)] = 5'^{x — y) or derivatives of the canonical 
fields may appear in their Lax operators violating the ultralocal condition and bringing additional 
complicities, which might not always be resolved. Due to this reason the theory and application for 
the nonultralocal models are still in the process of development and are far from completion. In 
spite of many important models belonging to this class, it is rather disappointing to note that, this 
category of models has not received the required attention in the literature. 

2 Integrable structures in ultralocal models 

We focus first on the ultralocal systems due to their relative simplicity and formulate a unifying 
scheme for generating such quantum and statistical integrable models. For ensuring the integrability 
of an ultralocal model it is sufficient to impose certain matrix commutation relation known as the 
quantum Yang-Baxter equation (QYBE) on its representative Lax operator in the form 

^o6(A - n)Laj{X)Lbj{n) = Lbj{n)Laj{X)Rab{>< - 1^), (2.1) 

defined at each lattice site j = 1,2, ... ,N. The above QYBE expresses actually the commutation 
relations among different matrix elements of the L-operator, given in a compact matrix form, where 
the structure constants are determined by the spectral parameter dependent c-number elements of 
the R{X — /x)-matrix. The i?-matrix in turn should satisfy a similar but simpler YBE 

RabiX - ll)Rac{\ - l)Rbcil^ - 7) = Rbci^i - 7)Kc(A - l)Rab{\ " M)- (2-2) 

Since our intention is to establish the integrability which is a global property, we have to switch from 
this local picture at each site j to a global one by defining a matrix, known as the monodromy matrix 

r„(A)^n.,,A). (2.3) 

Multiplying therefore the QYBE (2.1) for j = 1, 2, . . . , A'^ and using the ultralocality condition (1.2), 
thanks to which one can treat the objects at different lattice points as commuting objects as in the 
classical case and drag Laj{X) through all Lbfc(A*)'s for k ^ j,h^ ato arrive at the global QYBE 

Rab{X- ^x)Ta{X)n{^x)=n{^x) Ta{X) RabiX- fi). (2.4) 

Note that the local and the global QYBE have exactly the same structural form. Invariance of the 
algebraic form also for the tensor product of the algebras, as revealed here, indicates the occurrence of 
the coproduct related to a deep Hopf algebra structure underlying all integrable systems [4] . We will 
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see below that for nonultralocal models such a structure is modified a bit to include additional braiding 
relations. For the periodic ultralocal models defining further the transfer matrix as r(A) = traTa{X), 
taking trace from both sides of the global YBE (2.4) and canceling the i?-matrices due to the cyclic 
rotation of matrices under the trace we reach finally for r(A) at the trace identity (1.1) defining 
the quantum integrability of the system. Therefore we may conclude that the local QYBE (2.1) in 
association with the ultralocality condition (1.2) is the sufficient condition for quantum integrability 
of an ultralocal system. Consequently we may define such an integrable system by its representative 
Lax operator together with the associated i?-matrix satisfying these criteria. Note that we are 
concerned here only with the systems with periodic boundary condition. For models with open 
boundaries, the QYBE should however be modified with the inclusion of a reflection matrix, which 
was introduced in detail in [5]. 

2.1 List of well known ultralocal models 

To have a concrete picture before us, we furnish a list of the well known ultralocal models together 
with their L-operators and i?-matrices. We will however restrict here for simplicity only to the 
quantum models with 2 x 2-matrix Lax operators associated with 4x4 i2-matrices. We show in 
the next section how these Lax operators can be generated in a systematic way confirming their 
integrability. The i?"^-matrix that satisfies the YBE relation (2.2), with the indices taking values 

1. 2 only, can be given in a simple form by defining its nontrivial elements as [6] 

-^11 = -^22 = -^12 = -^21 = -^21 = -^12 = C. (2.5) 

These elements may be expressed explicitly through trigonometric functions in spectral parameters 
as 

a(A) = sin(A + a), 6(A)=sinA, c = sinQ; (2.6) 
or as its a — 0, A — limit, through rational functions as 

a(A) = A + a, 6(A) = A, c = a. (2.7) 
Moreover under a twisting transformation 

R{X) ^ R{X, 9) = F{e)R{X)F{e), with Fab{9) = e'^^^'-'^l^ (2.8) 



one gets twisted trigonometric and rational i?-matrix solution of (2.2) , which may be given by (2.5) 

with the difference = 6(A)e*^, i?2i = 6(A)e^*^. Apart from these /^-matrices there can be elliptic 

i?-matrix solution, for example that related to the XY Z spin chain and the 8- vertex model [1]. All 

models we consider here however are associated with the trigonometric or the rational ii-matrices 

and in the list presented below we group them accordingly, denoting H for the Hamiltonian and 

C {Ln) for the Lax operator related to field (lattice) models. 

I. Models associated with trigonometric i?-matrix (g = e'",^ = e*^) 

i) Field models 

1. Sine-Gordon model [7] 

TV? ^ ( Wi 'fn sin(A — au) \ 

Utt-Uxx = — sm(au), L=\ \ p = u. (2.9) 

OL \ msm(A + auj, —ip I 



4 



2. Liouville model [8] 

utt - u,. = e*°", c = t^ i^iaJ'^'Zp ) ' i<^)^piy)] = i^i^ - y)- (2-10) 

3. A derivative NLS (DNLS) model [9] 

^V'^-V'.. + 4^VVV'. = 0, C = i\ ), iV = VV, [^(x), = <5(x-y) (2.11) 

4. Massive Thirring (bosonic) model (MTM) [6] 



- < ^^^^^^^ ) ■ ^ - - ^^"'-^ - (-^) 

ii) Lattice Models 

1. Anisotropic spin chain [10] 

N 

H = Y^ + C7^cr^+i + cos ao-^cT^+i, Ln(6 = sin(A + aaV^) + sina (a+(7~ + (t~(t+) (2.13) 

n 

2. Lattice SG model [11] 

^ /,N / '/(itn) e'*'"^, mAsin(A — ait„) \ o, , 9,9 , ^ / 

^"*^'=(rAsi„(A + ™„), .-U «=(«») = l + ".^A^cos<.(2^ + l) (2.14) 

3. Lattice Liouville model [8] 

= ( VKfe-C ) . = ' + A=.'«(-.«>. (2.15) 

4. Lattice DNLS model [12] 

T iP\ Mr^" - ^At \ , cosa(2jV„ + l) 

^"(0 = ^^^^ + ^-(iv„+i) J ' K,^L] = ^^^^ . (2.16) 

5. Lattice MTM [13] 

Exact lattice version of MTM (2.12). 

Lax operator: L„ = Ln^Ln^ (each factor is a realization of (2.16) for a bosonic mode). 

6. Discrete-time or relativistic quantum Toda chain [14] 

H = Y^ (cosh2api + cosha(pi +pi+i)e("^-"'+i)) ^ = j f 



-ae""", 



(2.17) 
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la. Models associated with twisted trigonometric i?-matrix 

6. a) Quantum Suris discrete-time Toda chain [15, 14] 



-ae 



7. Ablowitz-Ladik model [16, 6] 



(2.18) 



ibj,t + (1 + ab]bj)ibj+, + = 0, LkiO = ^ 



, [bkM] = ^klil - bibk) (2.19) 



II. Models associated with rational i?-matrix 

i) Field models: 

1. Nonlinear Schrodinger equation (NLS) 

+ ^xx + (V'V)V' = 0, jO.{\) -- 

ii) Lattice Models: 

1. Isotropic XXX spin chain [10] 

^ ^ . 1 ^ ^ 

-H' = X]^" ■ "^"+1' -^an (A) = AI + QfPan, -Pan = -^(I + <?« " 



A, V \ 

Vt, -A i ■ 



(2.20) 



(2.21) 



2. Gaudin model [18] 

In the simplest case the Hamiltonians 



N 



1 



Hk = Y,^—{dk-di), k = l,2,...,N, Lak{X) = {X-ek)I + aPak. (2.22) 



l=^k 



3. Lattice NLS model [11] 

L(\)=( ^ + A-2{2s-ANn)hl 

4. Simple lattice NLS [19] 



,Nn = i^li^n, [A,i^l]=Skl- (2.23) 



5. Discrete self trapping dimer model [20] 



H 



i j^isa - iV(«))2 + (V't(l)V.(2) ^ ^t(2)^(l)) 



(2.24) 



,[^^^\^^^'^]=6ab,a,b=l,2 (2.25) 



Lax operator L(A) = L^^\X)L^'^\X), ( each factor as (2.24) for each of two bosonic modes). 
6. Toda chain (nonrelativistic) [6] 



^=y:(k+e(— )V MA) = ( ^"-^ 



I 







(2.26) 
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3 Unifying algebraic approach in ultralocal models 



Though the QYBE itself represents an unifying approach for all ultralocal models, we intend to 
specify here a common algebraic structure independent of the spectral parameter that will not only 
systematize the models including those listed above, but also identify their common integrable origin, 
establishing naturally the quantum integrability for all of them, simultaneously. From the above list 
of models, one may observe that, different integrable models have their representative Lax operators 
in diverse forms with varied dependence on the spectral parameter as well as on the basic operators 
like spin, bosonic or the canonical operators. However the i2-matrices associated with all of them are 
given by the same form (2.5) with known trigonometric (2.6) or its limiting rational (2.7) solutions. 
To explain this intriguing observation we may look for a common origin for the Lax operators linked 
with a general underlying algebra free from spectral parameters, though derivable from the QYBE. 
We propose to take the Lax operator of such ancestor model in the form [17] 

4.. (0=^ ^4e-^-SS^^-i^^.aSS )^ ^ = e ,e± = 2sma^ , (3.1) 

where 5^ and c^, a = 1. 2 arc some operators, the algebraic properties of which arc specified below. 
The structure of (3.1) becomes clearer if we notice the decomposition L['l!^g\^) = ^L^ + ^^^L_, 
where L± are spectral parameter ^-independent upper and lower triangular matrices similar to the 
construction of [4]. Inserting (3.1) in QYBE together with its associated i?-matrix (2.5) with trigono- 
metric solution (2.6) and matching different powers of the ( we obtain the underlying general algebra 
as 

[5^ 5^] = ±S^, [S+, S~] = (m+ sm{2aS^) + M~ cos(2a5^)) -J—, [M^ , ■] = 0, (3.2) 

with = ib^-v/±T(c^C2 ± c^c^) behaving as central elements with arbitrary values of c's. As 
we have mentioned above, the integrable systems are associated with an important Hopf algebra A, 
exhibiting the properties like 1) coproduct A{x) : A ^ A A, 2) antipode or 'inverse' S : A ^ A, 3) 
counit e : A ^ k, A) multiplication M : Ai^ A ^ A and 5) unit a : k ^ A. It can be shown that all 
these properties hold also for (3.2) defining it as a Hopf algebra. Referring the interested readers to 
the original works [21] for more mathematical treatment of the noncocommutative Hopf algebra, we 
give here only some simple and intuitive arguments in its constructions. For example the coproduct 
A(x), the most important of these characteristics, can be derived for algebra (3.2) by exploiting a 
QYBE property that the product of two Lax operators LajLaj+i is again a solution of the QYBE 
and may be given in the explicit form as 



A(5+) = c^-e'"^ (8) S+ + 5+ ® cje-^"^ , A(5-) = c^e'"^ (g) 5" + (g) 



e 



A{S^) = I®S^ + S^®I , A{cf) = cf^cf. (3.3) 

The multiplication property mentioned above is also in agreement with the ultralocality condition, 
which is used for transition from local to global QYBE following the multiplication like 

{A^B){C^D) = {AC^ BD) with A = Lj (A) , 5 = (/x) , C = L^+i (A) , L» = L^+i (/x) . (3.4) 
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Note that (3.2) is a g-deformed algebra and a generalization of the well known quantum algebra [21] 
Ug{su{2)). 

In fact different choices of the central elements reduce this algebra to the q-spin, q-boson 
as well as various other q-deformed algebras along with their undeformed limits. Therefore we can 
obtain easily the coproduct for these algebras, whenever admissible, from their general form (3.3) in 
a systematic way by taking the corresponding values of d s. 

3.1 Generation of models 

We know that the well known integrable models listed above were discovered at different points of 
time, mostly in an isolated way and generally by quantization of the existing classical models. How- 
ever, as we will sec, they can actually be generated in a systematic way through various realizations 
of the same Lax operator (3.1) giving a unifying picture of integrable ultralocal models. For this we 
find first a representation of (3.2) like 

= u, S+ = e-'Pg{u), 5" = g{u)e'P, (3.5) 

in physical variables with [u,p] = i, where the operator function 

g{u) = {k + sma{s - u){M'^ sma{u + s + 1) + M~ cosa{u + s + 1)))^ -r^, (3.6) 

Sin ex 

containing free parameters k and s. We demonstrate now that the Lax operator (3.1), which rep- 
resents a generalized lattice SG like model for (3.5) may serve as an ancestor model (with possible 
realizations also in other physical variables like bosonic ip,il^^ or spin s='=,s^ operators) for gener- 
ating all integrable ultralocal quantum as well as statistical systems. As an added advantage, the 
Lax operators of these models are derived automatically from (3.1), while the i?-matrix is simply 
inherited. The underlying algebras of the models are also given by the corresponding representations 
of the ancestor algebra (3.2), which being a direct consquence of the QYBE ensures the quantum 
integrability of all its descendant models, that we construct here. It should be stressed that due to 
the symmetry of the solution (2.5): [R{X — IJ^),cr"' fS> cr"'] =0, a = 1,2,3 the Lax operator (3.1) as 
a solution of QYBE may be right or left multiplied by any a". We shall use this freedom in our 
following constructions, whenever needed. 

Note that we may generate also the quantum field models by taking properly the continuum limit 
of their lattice variants with the lattice spacing A — 0. Though in general such transitions to the 
field limit might be tricky and problematic we suppose their validity assuming the lattice operators 
to go smoothly to the field operators pj — >■ p{x),jpj ipix), with the corresponding commutators: 
[il^j, t/jj,] = ^Sjk ii^i^), i^{y)] = S{x — y) etc. The lattice Lax operator therefore should reduce to 
its field counterpart C{x,\) as Lj{X) I -\- iAC{x,X) + O(A^). The associated i?-matrix however 
remains the same, since it docs not contain lattice constant A. Thus integrable field models like 
sine-Gordon, Liouville, NLS or the derivative NLS models can be recovered from their exact lattice 
versions and having the same quantum i?-matrix, though all discrete models may not always have 
such a direct field limit. 
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3.1.1 Models belonging to trigonometric class 

1.) Choosing trivially all central elements as = 1, a = 1,2, which gives M~ = 0, M+ = 1, (3.2) 
reduces clearly to the well known quantum algebra Uq{su{2)) [21] given by 

[5^5±]=±5±, = [253]g. (3.7) 

with the known form of its coproduct recovered easily from (3.3). The simplest representation 
S = ■^a for this case derives from (3.1) the integrable XXZ spin chain (2.13). On the other hand, 
representation (3.5) with the corresponding reduction of (3.6) as g{u) = jiHTo +cosa{2u+ 1)]^ 
with suitable choice of parameters s, k recovers the Lax operator of lattice sine-Gordon model (2.14) 
directly from (3.1) and at its field limit the field Lax operator (2.9). Note that the spectral dependence 
in e± appearing in (3.1) can be easily removed through a simple gauge transformation [4] and therefore 
we ignore them in our construction and use the freedom of translational symmetry of the spectral 
parameters A — ^ A + const., whenever needed. 

2. ) An unusual exponentially deformed algebra can be generated from (3.2) by fixing the elements 
as cf = C2 = 1, Ci = ^2 = 0, which gives = ±^^/±l and 

[^'•**i=±«' l«"-»"l = 2?^ (3.8) 

and reduces (3.6) to g{u) = 72sma " '^^^^ algebra and the corresponding realization yields 
clearly from (3.1) the Lax operator of the lattice Liouville model (2.15) and at its field limit that of 
the Liouville field model (2.10). 

It is interesting to observe here that though the underlying algebraic structure and hence its 
realization are fixed by the choice of M^, the Lax operator (3.1) depends explicitly on the set of c's 
and therefore may take different forms for the same model. For example in the present case ^ 
would give again the same value for but a different Liouville Lax operator [22] more convenient 
for the Bethe ansatz solution. 

This opens up therefore interesting possibilities for obtaining systematically different useful Lax 
operators for the same integrable model, as well as for constructing new nonultralocal models [23]. 

3. ) Recall that the well known q'-bosonic algebra may be given by [24] [^,iV] = A, [A^,N] = 
—A^, AA^ — q~^A^A = q^^ or in its conjugate form with q — ^ q~^. Combining these two forms we 
can easily write the commutator of such q-bosons as 



[A,N]=A, [A\N] = -A\ [A^t] = ^»K"(2^ + l)). 



cos a 



(3.9) 



It is interesting to find that for the choice of the central elements cf = C2 = 1, c^ = —iq, c^ = | 
compatible with M+ = 2sin a, M" = 2icos a we may get a realization 

S+ = -kA, S' = kA^, = -N, K = -z(cotQ)^, (3.10) 

with (3.2) reducing directly to the relation (3.9), which gives thus a new integrable q-boson model. It 
is important to note now that either using (3.5) which simplifies (3.6) to g'^{u) = [—2u]q or directly 
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taking the mapping of the (7-bosons to standard bosons: A = V'( 2^0050 )^ ' ^ ~ "^"^"^ ' may 
convert (3.1) with (3.10) to an exact lattice version of the quantum derivative nonlinear Schrddinger 
(QDNLS) equation (2.16) and consequently to the QDNLS field model (2.11). The QDNLS is related 
also to the interacting bose gas with derivative 5-function potential [25]. 

4. ) Since the matrix product of Lax operators with each factors representing different Lax operator 
realization for the same model should give again a QYBE solution, we can construct multi-mode 
integrable extensions by taking the product of single-mode Lax operators. Using this trick, i.e. by 
combining two QDNLS models constructed above as L{c^ ,ip^^^)L{c2 , cf ,tp^'^^) = L{X), we can 
create further an integrable exact lattice version of the massive Thirring model [13]. At the continuum 
limit it goes to the bosonic massive Thirring model introduced in [6], the field Lax operator (2.12) of 
which can be given simply by the superposition jC = £^^\^,k±,ip^^^) + a^£^'^\^,kzf,ip^'^^)a^, where 
1 lb ik± sina = e^'a' and the constituing operators jC^^") is given clearly by the DNLS Lax operator 
(2.11) for each of its two bosonic modes. 

5. ) Since the general algebra permits trivial eigenvalues for central elements, one may choose both 
= 0, which might correspond to different sets of choices like i) = 1 , a = 1, 2, or ii) = 

1 , a = 1, 2, or Hi) = ±1, or iv) = 1, with the rest of c/s being zero. It is easy to see that 
all of these sets lead to the same underlying algebra 

[S+, S-] = 0, [5^ S^] = ±S^, (3.11) 

though generating different Lax operators from (3.1). 

As here (3.6) gives simply g{u) =const., interchanging canonically u — > —ip,p — —iu, from (3.5) 
one gets 

= -ip, = ae^" , (3.12) 

which evidently generates from the same general Lax operator (3.1) the discrete-time or relativistic 
quantum Toda chain (2.17). Note that, iii) and iv) give two different Lax operators found in [14] and 
[26] for the relativistic Toda chain. Case i) and ii) on the other hand could be used for constructing 
nonultralocal quantum models, namely light-cone SG and the mKdV model [23]. 

Models in twisted trigonometric class: Under twisting when the i?-matrix changes as (2.8) 
the associated Lax operator is also transformed similarly as -Ln(A) Lni\,0) = Fn{6)Ln{X)Fn{6), 
with Fan{0) = e^^^'^a-Sn) _ As a result the ancestor model (3.1) associated with the trigonometric 
twisted i?-matrix (2.8) gets deformed with its operator elements changing as 

4-c^e-^"s ^ Si = e-^^^'t Si (3.13) 

and as a consequence the diagonal elements of the twisted Lax operator take the form e*(^='=")'^fc , with 
obvious preferance for the choice 9 = ita. 

6. ) We may generate the quantum analog of Suris discrete-time Toda chain belonging to the 
twisted class by starting from the ancestor model with the change (3.13), but by fixing the parameter 
9 = —a (an equivalent model is obtained by the choice 9 = a). Using the same realization (3.12) we 
arrive now at the explicit form (2.18). 

7. ) However if we start from the same twisted ancestor model with the same value 9 = —a of 
the twisting parameter, but take the central elements as = = with c]" = C2" = 1 giving 
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= ^\/±T (compare with the Liouvillc case!), all noncommuting operators clearly vanish from 
the diagonal elements of the resulting Lax operator. Moreover renaming the deformed operators as 
bk = 2sina5^, we get their modified algebra as a type of g-boson : [&fc,b]^] = 5^/(1 — b^bk) and 
thus generate finally the exact form of the Ablowitz-Ladik model (2.19). 

The domain of the models considered can therefore be considerably extended if we use twisting 
and some other allowed transformations [27] that preserves integrability. 

3.1.2 Models belonging to the rational class 

One of the crucial parameters inbuilt in both i?-matrix and above Lax operators is the deformation 
parameter q = e*", the physical meaning of which is as anisotropic or relativistic parameter. We 
consider now the undeformed limit q ^ 1 or a ^ related to isotropic or nonrelativistic models 
belonging to the rational class, which reduces various a-dependent objects as is^,{c^} 
{c^}, M"*" —m^,M~ —am^, ^ ^ 1 + iaX. This transforms (3.2) to a ^-independent algebra 
with 

[s+, s~] = 2m+s^ + m-, [s^, s±] = (3.14) 

were m"'" = c5c2, m~ = c\c2 + CiC^ and c^,i = 0,1 are central to (3.14). Note that (3.14) is a 
generalization of spin as well as the bosonic algebra and its coproduct can be obtained as a limit of 
(3.3). Consequently, the general Lax operator (3.1) is converted into 

and the quantum i?-matrix (2.5) is reduced to its rational form (2.7). 

We would see that the ultralocal integrable systems belonging to the rational class can be gen- 
erated in a similar way now from the Lax operator (3.15) with algebra (3.14), all sharing the same 
rational i?-matrix (2.7). It is not difficult to check by a variable change {u,p) that at the 

limit a — s- (3.5) reduces to a generalized Holstein-Primakov transformation (HPT) 

s^ = s-N, s+=go{N)iP, s- = il;^go{N), gl{N) = m- + m+{2s - N), AT = y^ty,^ (3J6) 

which is also an exact realization of (3.14). Therefore Lax operator (3.15) with such a realization may 
be considered as a generalized lattice NLS, which would serve as a generating model for all quantum 
integrable models belonging to the rational class. 

1.) The choice m+ = 1, m~ = 0, clearly reduces (3.14) to su{2) algebra [s"*", s"] = 2s^, [s^, s^] = 
±s^. A compatible choice = l,c^ = yields from (3.15) for the spin-^ representation the Lax 
operator of the XXX spin chain (2.21). 

Taking spin-^ and spin-1 realizations alternatively along the lattice we can construct now the 
integrable alternate spin model discovered in [29]. 

Note that a slightly different choice = — C2 = 1, = giving m+ = — 1, m~ = 0, generates on 
the other hand the corresponding model with su{l, 1) algebra. 

The bosonic realization (3.16) in present cases with m'^ = ±l,m~ = 0, is simplified to the 
standard HPT with 5o(-^) — i(2s — V'V)) which reproduces from (3.15) the exact lattice NLS 
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model (2.23) and at the continuum limit the more familiar NLS field model (2.20), with +(— ) sign 
in the HPT corresponding to the attractive (repulsive) interaction. 

2. ) A complementary choice m"'" = 0, m~ = 1, on the other hand converts (3.16) to = 
ijj, s~ = , = s — N due to goiN) = 1 and reduce (3.14) directly to the standard bosonic relations 
[■0, AT] = [ijj^,N] = —il^^, [V', V'^] = 1- Remarkably, (3.15) with this realization generates yet 
another simple lattice NLS model with Lax operator (2.24). 

3. ) Combining two such bosonic Lax operators (2.24), constructed above: L(^)(A)L(^)(A) = L{X) 
and considering them to be inserted at a single site we can construct the Lax operator of an integrable 
model involving two-bosonic modes, which yields the quantum discrete self trapping model (2.25). 

4. ) Note that the trivial choice = gives again algebra (3.11) and hence the realization 
(3.12). This however yields from (3.15) the Lax operator of the nonrelativistic Toda chain (2.26) 
associated with the rational i?-matrix. It is interesting to note that in [30] the Lax operators like 
(3.1) and (3.15) appeared in their bosonic realization and were shown to be the most general possible 
form within their respective class. 

Therefore these quantum Lax operators are in the core of the ultralocal integrable models, both 
discrete and continuum, which can be constructed from them in a unified way. Models belonging to 
trigonometric and rational class arc generated from (3.1) and its limiting form (3.15) respectively, 
and therefore inherit the same corresponding i?-matrices (2.6) and (2.7) . 

3.2 Fundamental and regular models 

The Lax operator Laj(X) in general acts on the product space Va ® hj, of the common auxiliary 
space Va and the quantum space hj at site j. The models with Va isomorphic to all hj,j = I, ... N 
and given by the fundamental representation are called fundamental models. For such models the 
finite-dimensional matrix representations of auxiliary and quantum spaces become equivalent and 
may lead to Lai {^) = Ral (-^) • However for clarifying a misconception prevelant in the literature, we 
should stress that a model is represented by its Lax operator only, the associated i?-matrix accounts 
for the commutation property of the elements of this Lax operator through QYBE. Therefore even 
for a fundamental model the Lax operator may differ from its i?-matrix. We may however demand 
for the fundamental models an useful additional property: Lai{0) = Pal, known as the regularity 
condition given through the permutation operator, which may be expressed in the general case as a 
X matrix P^"-* = J2^a=i ^afs-^jda^ being a matrix with its (a,/?) element as 1 and the rest 

0. 

Recall that the global operator r(A) is constructed from the local Lax operators Laj{\), j = 
1, . . . N as r(A) = tra {Lai{\) ■ ■ ■ LaNW) , where the transfer matrix t(A) acts on the total quantum 
space Ti = Therefore, from the knowledge of Lax operators it is possible to derive all 

conserved quantities including the Hamiltonian, which in general would be nonlocal objects. The 
above regularity condition on Lax operators however allows to overcome this difficulty and obtain 
Hamiltonians with nearest-neibour (NN) interactions. Let us pay some special attention to this 
specific group of models since, as we will see below, the most important integrable models applicable 
to condensed matter physics related problems are given by the regular models with n = 2,3,4 etc. 
For this reason, though our main concern in this paper is 2 x 2 auxiliary matrix space, we describe 
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here more general n cases and demonstrate that Hamiltonians of such different physical models 
interestingly have similar forms, when expressed through the permutation operator Pjj+i- Such a 
permutation operator exhibits space interchanging property PajL^k = Lj^Paj, along with P^ = 1 
and tra{Paj) = 1- For all regular and periodic models, using the freedom of cyclic rotation of matrices 
under the trace, we can express the transfer matrix as 

r(0) = tra (PajPaj+l.-PaNPal-Paj-l) = {Pjj+L.PjN Pjl..Pjj-l)traiPaj) (3.17) 

for any j and as its derivative with respect to A we similarly get 

N N 
r\0) = traY.(PajL'^j+M-PaNPal..Paj-l) = Y.i^j+li^)--PjNPjl..Pjj-l)tra{Paj), (3.18) 

where the periodic boundary condition: LaN+j = Laj is assumed. Defining now H = ci = 
^lnT(A)|;^=o = and using (3.17), (3.18) we may construct the related Hamiltonian 

as 

N 

^ = E L'jj+ii^) Pjj+i (3-19) 

with only NN interactions, where all nonlocal factors are canceled out due to relevant properties of 
the permutation operator. Similarly taking higher derivatives of Inr(A), higher conserved quantities 
Cj,j = 2,3, ...N can be constructed for these regular models. Note that the conserved operator cj 
involves interactions of j + 1 neighbors. 

For the simplest case of n = 2, we may take the Lax operator as the i?-matrix given by (2.5), which 
satisfies clearly the regularity condition L{0) = R{0) = P for both trigonometric and rational cases. 
Moreover for (2.6) the part L'jjj^^{0) in (3.19) introduces anisotropy reproducing the Hamiltonian 
of the XXZ spin chain (2.13). However since for the rational case (2.7) L'jj^^{0) = 1, using the 
expression Pjf+i = E|a=i ^ip^jkx^ = + where -ffjj+i = ^j^j+i, (3.19) is reduced clearly 

to the isotropic spin-i Hamiltonian H"^ = J2j ^jj+i (^-^l)- 

It is intriguing to note that for the rational class the same form of Hamiltonian H = J^jLi Pjjli 
with several higher values of n describes most of the important integrable models, though their 
physical forms are given mainly through various representations of the permutation operator. For 
example, for n = 3 corresponding to SU{3) group we can express the permutation operator Pjj\.i = 
J2^a=i ^a/s^jsti^ through spin-1 operators S giving a variant of the integrable spin-1 model 

H = J2SjSj+i + e{SjSj+if, withe = +l. (3.20) 

3 

Considering a supersymmetric invariant ^^(1,2) case, i.e. realizing the corresponding graded permu- 
tation operator Pjj'^l using fermionic {caj,clj), a =|, J, and spin S operators we may construct again 
from the Hamiltonian density Hjj_^_i = {2Pjj!^l — 1) the well known integrable t — J model [31] 



H 



tJ 



= E ^JJ+i = E f E + h.c. ]P+ J (SjSj+i - jnjnj+i) + nj + n,+i (3.21) 
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with J = 2t = 2, where V projects out the double occupancy states. 

A different 4-dimensional realization of the fermion operators on the other hand converts the 
same Hamiltonian to an integrable correlated electron model proposed in [32]. 

Similarly for n = 4, i.e. for SU{A), realizing P jf j^i = Pjf+i ® ^jf+i factorized form we can 

get 

^=E^;li®^;ii = iE(^j^+i+i)(^L+i+i) = i[^'^+^^+E(^^+i^L+i+i)] . (3-22) 

3 j j 

with H'^''^ representing isotropic spin-^ Hamiltonians (2.21). Adding now interaction along the rung: 
Hrung = ^j^ji with [H^Hmng] = to (3.22), where £7, r represent the spins along two legs of 
the ladder, we may construct a model which is nothing but the integrable spin-i ladder discovered 
recently [33]. 

On the other hand, from the same form of Hamiltonian but by considering a supersymmetric ex- 
tension SU{2, 2) we may realize -Pjl+i again through fermion operators (caj, c^^), a =t, j, to construct 
an integrable extension of the Hubbard model proposed in [34]. 

One can repeat the above construction of the spin-ladder model for generating also an integrable 
t — J ladder model introduced in [35], which would therefore corresponds to a similar construction 
in n = 6 with Hamiltonian 

3 3 3 

= l[if(i)*^ + H('^^^ + + 1)]' (3-23) 

3 

where i?^")*"^,a =1,2 arc t — J Hamiltonians (3.21) along two legs. Adding a suitable Hmng to 
(3.23) with [H, Hmng] = 0, defining the interaction along the rung we finally obtain the integrable 
t — J ladder model. 

Apart from the above applications of integrable systems having similar structure from the alge- 
braic point of view, we should mention some other important models like Hubbard model and the 
Kondo problem, which also falls in the class of exactly solvable problems in one-dimension [36] . Em- 
ploying further twisting and gauge transformations on multi-fermion or multi-spin integrable models 
one can generate another type of integrable models of current interest [37]. Importance of solvable 
models in physical systems, their relevance to experiments and related isuues are discussed in [38]. 
For detailed and involved application of Bethc ansatz technique including that for the theory of cor- 
relation functions to various integrable systems like 5-bose gas, NLS, sine-Gordon etc. the readers 
are referred to [3]. 

3.3 Fusion method 

We have constructed spin, boson as well as the q-spin and q-boson models through realization of 
particular Lax operators with inequivalent auxiliary and quantum spaces. However in case of finite- 
dimensional higher rank spin representations there exists an intriguing method, known as the fusion 
method, for obtaining higher spin models by fusing the elementary i?-matrices like (2.5). Thus by 
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fusion of only the quantum spaces one can construct spin-s Lax operators with a spin-| auxiliary 
space, as obtainable also directly from (3.1) as a particular realization. Fusing further the auxiliary 
spaces the higher-spin Lax operator with spin-s auxiliary space may be constructed as 

Lab = {P^ ® P^) n n ^aM^ + -f"- ■^■))(^a+ ® Pb)92sW (3.24) 

j=l k=l 

with -f'a'^b) as the symmetrizer in the fused spin-s space a(b) and g2s some normalizing factor [39] . For 
the rational i?-matrix corresponding to (2.7) one obtains from (3.24) the intcgrable spin-s Babujian- 
Takhtajan model [39], which for s = I may be given in the same form as Hamiltonian (3.20), but with 
e = — 1. Similarly for the trigonometric case (2.6) the fused model would correspond to integrable 
anisotropic higher-spin chain. 

It may however be stressed that such fusion technique, as far as we know, has not been formulated 
yet for bosonic and q-bosonic models. Such extension , at least for the restricted values of q, needs 
therefore more attention. 

3.4 Construction of classical models 

The systematic procedure for constructing quantum integrable models as various reductions of the 

same ancestor model, as described here, is applicable naturally also to the corresponding classical 
models by taking the classical limit ?i ^ 0. At this limit all field operators would be transformed 
to ordinary functions with their commutators reducing to the Poisson brackets. Note also that 
parameter a appearing in the i?-matrix is scaled actually as ha, which yields the classical r-matrix: 
R{X) = I + hr{X) + 0{h'^) and reduces QYBE (2.1) to its classical limit {Lai{X),Lbj{iJ,)} = Sij[rabiX- 
/j.), Lai{X)Li,j{fx)]. The classical Lax operator reduced from (3.1) would remain however almost in 
the same form, though the corresponding quantum algebras would change into their corresponding 
Poisson algebras. This aspect of classical integrable systems is given in great detail in the excellent 
monograph [40]. Using these classical analogs of quantum systems one can therefore apply the 
algebraic scheme formulated above for generating quantum models also in classical context and 
construct systematically the corresponding classical integrable models [41]. 

4 Integrable statistical systems: vertex models 

D-dimensional quantum systems are known to be related to (1 -|- D)-dimensional classical statistical 
models, which is true naturally also for D = 1, where the integrability of models might get manifested. 
Interestingly, integrable quantum spin chain and the corresponding vertex model share the same 
quantum i?-matrix and have the same representation for the transfer matrix, commutativity of 
which: [t{X),t{h)] = guarantees their integrability. However, while the spin chain Hamiltonian Hg 
is expressed through the transfer matrix as Inr^X) = I + XHs + O(A^), the partition function Z of 
the vertex model is constructed from r(A) as Z = tr{T{X)^). The known integrable vertex models 
are usually related to the quantum fundamental models described above. 

In conventional vertex models each bond connecting N x M arrays in a 2-dimensional lattice can 
take n different possible random values with certain probabilities, which for a configuration i,j;k,l 
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of bonds meeting at each vertex point is given by the Boltzmann weights Wij^^i. These Boltzmann 
weights may be assigned as matrix elements Wij^klW = ^kiW of ^ -R-matrix (though it might be 
of a more general L-operator, as we will see below), which for integrable models must satisfy the 
Yang-Baxter equation (2.1) and correspond to a quantum integrable model. The partition function 
of these vertex models may be expressed as Z = J2configIla,b,j,k^a,j;b,kW- 

The simplest among the vertex models for n = 2 is the 6-vertex model [1], which corresponds 
to the XX Z spin chain and may be defined on a square lattice with a random direction on each 
bond ( left or right on the horizontal, up or down on the vertical), constrained by the ice rule, 
that the number of incoming and outgoing arrows at each vertex are the same. This leaves only 
6 possible configurations and the corresponding Boltzmann weights may be given by 6 nontrivial 
matrix elements of the J?-matrix (2.5) with (2.6). It is fascinating that this model may describe the 
possible configurations of Hydrogen (H) ions around Oxygen (O) atoms in an ice crystal having two 
different (close-removed) positions of the H-ions relative to the 0-atom in the H-bonding, while the 
ice rule corresponds to the charge neutrality of the water molecule. 

A more general 6-vertex model may be obtained if instead we assign its Boltzmann weights 
directly to the spin-^ matrix representation of the general ancestor Lax operator (3.1). The param- 
eters cf = —c^ = = — C2 = P- present in the Lax operator may be combined to serve as the 
horizontal h ~ lnp_|_p_ and vertical ~ In ^ fields acting on the model, which recovers amazingly 
the most general 6-vertex model proposed many years ago [42] through a different construction. This 
also confirms the fact that the Lax operator (3.1) is indeed in the core of integrable quantum as well 
as statistical models. Using twisting transformation one can recover also the 6^(1) vertex model 
introduced in [27]. 

We may consider higher vertex models with n > 2, which may be obtained from the i?-matrix 
(or the Lax operators) of the corresponding quantum integrable fundamental models with higher- 
dimensional auxiliary spaces. The well known examples are the 19-vertex model [43] related to the 
Babujian-Takhtajan integrable spin-1 model [39], the Boltzmann weights of which may be given by 
the matrix elements of Lax operator (3.24) with s = 1. Similarly one may construct the vertex 
models equivalent to the Hubbard model, supersymmetric t-J model, Bariev chain etc. [44]. 

In a following section new type of vertex models will be constructed from our ancestor Lax 
operator using nonfundamental representations. 

5 Directions for constructing new classes of ultralocal models 

The same unified scheme described in sect. 3 for constructing integrable models may be used also to 
indicate various directions for generating new integrable classes of quantum and statistical models. 

5.1 Inhomogeneous models 

In all above constructions the central elements in the ancestor models (3.1) or (3.15) are chosen as 
constant parameters. However if they are chosen as site dependent (or may even time dependent) 
functions we can get an inhomogeneous class of models. In these cases the c's would be attached with 
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site indices as c^s in the Lax operators and similarly in (3.6) Mj would appear as functions, leading to 
the corresponding inhomogeneous extensions of the known integrable models, namely inhomogeneous 
lattice sine-Gordon, Liouville, Toda chain, NLS model etc. However since the local algebra remains 
same as the original model, they have the same quantum i?-matrices. Though similar inhomogeneous 
Toda chain, NLS models etc. were proposed earlier as classical systems, they seem to be new and yet 
unstudied as quantum models. Recall that the impurity models proposed earlier [45] fall into this 
class and are obtained by a particular choice of inhomogeneous c^'s which amounts to a shifting of 
the spectral parameter. Implementing the same idea to the XXX spin chain we notice that, if in its 
constructing along with = 1 we choose C2 = —c\ = ej resulting again = l,m~ = 0, we get the 
same form of the Lax operator, but with a shift Lj(X — ej), resulting that of the Gaudin model (2.22). 
Similarly higher spin representations as well as su{l, 1) variant would yield other generalizations of 
the same model. The commuting set of Hamiltonians for the Gaudin model may be generated from its 
transfer matrix at the limit a — ^ [18] as Hj = — ppjyi r(A — ^ e,), j = 1,2, . . . ,N. Remarkably, 

" (life (^-'^k)) 

the Gaudin model may be mapped into the integrable BSC model, which is of immense contemporary 
interest [28]. 

Physically such inhomogeneities may be interpreted as impurities, varying external fields, incom- 
mensuration etc. 

5.2 Hybrid models 

Another way of constructing new models is to use different realizations of algebras (3.2) or (3.14) 
at different lattice sites, depending on the type of the i?-matrix. For example one may consider 
spin-i and spin-1 representations of su{2) at alternate lattice sites, which was realized actually in 
[29]. However we can build more general inhomogeneous integrable models by considering different 
underlying algebras and different Lax operators at differing sites. The basic idea is that the Lax 
operators representing different models that are descended from the same ancestor model and share 
the same i?-matrix can be combined together to build various hybrid models preserving quantum 
integrability. For example, we may consider fermion-boson or spin-boson interacting models by 
inserting alternatively spin-^ and bosonic (or q-bosonic) Lax operators at alternate sites. One of 
such physical constructions would be the celebrated Jayncs-Cummings model. It is possible also to 
construct some exotic hybrid integrable models, an example of which could be a hybrid sine-Gordon- 
Liouville model, where for x > it would follow the sine-Gordon dynamics, while for a; < the 
Liouville dynamics! 

5.3 Nonfundamental statistical models 

Vertex models, as mentioned, are described generally by the i?-matrix of a regular quantum integrable 
model. However one can construct a new class of integrable vertex models by exploiting a richer 
variety of nonfundamental systems, where we define the Boltzmann weights as matrix elements of the 
generalized Lax operator (3.1): L^J^(u) = u!a,j-b,k{u), with the use of the explicit matrix representation 
for the basic operators S^,S^ as 

< s,m\S^\m,s >= mSm,fh, < s,fh\S'^\rn, s >= f^{m)Srn±i,fh- (5.1) 
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Here ff{m) = fj{m + 1) = g{m) is defined as in (3.6). Such general Boltzmann weights would 
now represent an ancestor vertex model analogous to the quantum case and would generate through 
various reductions new series of vertex models, linked to q-spin and q-boson with generic q, q roots 
of unity and g ^ 1 [46]. In all these models, generalizing the usual approach the horizontal (h) 
and vertical (v) links may become inequivalent and independent at every vertex point. The h links, 
which are related to the auxiliary space admit 2 values, while the v links, which correspond to 
the quantum space may have richer possibilities with j, fc G [li-D]) D being the dimension of the 
nonfundamental matrix-representation of the q-algebras. The familiar icc-rulc is generalized here 
as the 'color' conservation a + j = b + k for determining nonzero Boltzmann weights. Note that, 
alternatively finite-dimensional higher spin and g-spin vertex models can also be constructed using 
the fusion technique [39]. 

An interesting possibility of regulating dimension for the matrix representation opens up at = 
±1, when a variety of new q-spin and q-boson vertex models with finite-dimensional representation 
can be generated [46] . 

As in quantum models we can also construct here a rich collection of hybrid models by combining 
different vertex models of the same class and inserting their defining Boltzmann weights along the 
vertex points / = 1,2, ... in each row, in any but in the same manner. Due to the association 
with the same i2-matrix the integrability of such statistical models is naturally preserved. 

6 Unified Bethe ansatz solution 

In physical models our aim usually is to solve the eigenvalue problem for the Hamiltonian only. 

Solvable models allow such exact solutions H \ m >= Em \ m > through coordinate formulation of 
the Bethe ansatz (CBA) [47], which was used successfully in many condensed matter physics related 
problems like spin-chain, attractive and repulsive (5-Bose gas, Hubbard model etc [48]. Nevertheless 
CBA depends heavily on the structure of the Hamiltonian of individual models and lacks consequently 
the unified approach of its algebraic formulation. We would focus here briefly only on the algebraic 
Bethe ansatz (ABA) [2, 3], which under certain conditions can solve the eigenvalue problem for 
the spectral parameter-dependent transfer matrix t(A) | m >= A^(A) | m > and hence through its 
expansion the eigenvalue problem for the whole set of conserved operators, simultaneously. Moreover, 
the ABA due to its predominantly model-independent features, which we will demonstrate below, 
appears to be a fairly universal method. 

Since the eigenvectors are common for all commuting conserved operators, by expanding In Aot(A) 
simply as 

ci I m >= A^(0)A;;^(0) I m >, C2 | m >= (A'„(0)A„^(0))' | m > (6.1) 

etc. we obtain their respective values, where one may take H = ci or other combinations of c's as the 
Hamiltonian, depending on the concrete model. This powerful method applicable to both integrable 
quantum and statistical systems requires however explicit knowledge of the associated Lax operator 
and the i?-matrix. 

It may be noticed that the Oj^^-diagonal element B{X) (C(A)) of the monodromy matrix (2.3) acts 
generally like creation (annihilation) operator for the pseudoparticles, induced by the local creation 
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(annihilation) operator as the matrix elements in Lj{X) acting on the quantum space at j. Therefore 
the m-particle state | m > may be created by acting m times with B{Xa) on the pseudovacuum 
I >= nf I >j , giving I m >= B{Xi)B{X2) ■ ■ ■ B{Xm) \ >, where we suppose the crucial 
annihilation condition C(Aa) | >= 0. 

Now for solving the eigenvalue problem of t(A) = ^(A) + D{X) exactly, we have to drag this 
operator through the string of -B(Aa)'s without spoiling their structure and finally hit the pseudovac- 
uum giving A{X) \ >= a{X) \ > and D{X) \ >= /3(A) | >. For this purpose therefore one 
requires commutation relations between the elements of (2.3), which for ultralocal models may be 
derived from the QYBE (2.4). This, apart from ensuring the integrability of the system, is another 
important role played by (2.4), yielding the relations 

A{X)B{Xa) = f{Xa-X)B{Xa)A{X)-hiXa-X)B{X)A{Xa), 

D{X)BiXa) = f{X-Xa)B{Xa)DiX)-h{X-Xa)BiX)D{Xa), (6.2) 

together with the trivial commutators [A{X),A{n)] = [B{X),B{n)] = [D{X),D{n)] = [A{X),D{iJ,)] = 
etc., where /(A) = /i(A) = are combinations of the elements from the i?(A)-matrix (2.5). 
We notice that (6.2) are almost the right kind of relations but for the second terms in both the RHS, 
where the argument of B has changed spoiling the structure of the eigenvector. However, if we put 
the sum of all such unwanted terms = 0, we should be able to achieve our goal. In field models such 
unwanted terms vanish automatically, while in lattice models their removal amounts to the Bethe 
equations, which may be induced independently by the periodic boundary condition, giving 

"(M^^ rr /(Aa-Afe) 

Therefore the ABA finally solves the eigenvalue problem for t(A) yielding 

(rn \ / m \ 

n /(Aa - A) j «(A) + /(^ - j /^W' (6-4) 

where the Bethe equation (6.3), which is equivalent also to the singularity- free condition of the 
eigenvalue (6.4) serves in turn as the set of equations for determining the parameters A^. 

Note that in both the above equations a(A) = (< 0|L]i(A)|0 >)^ and /3(A) = (< 0|Lf (A)|0 >)^ 
are the only model dependent parts given by the action of the upper and lower diagonal operator 
elements L"(A),i = 1,2 of the Lax operator of the model on the pseudovacuum. For vertex models, 
for which the ABA formulation goes parallelly, the Lax operator elements in the above equations 
should be replaced by their matrix representations expressed through the Boltzmann weights as 
< 0|L]^(A)|0 >= cj+,i;+,i(A), < 0|L|^(A)|0 >= u;_^i;_^i(A). It is remarkable that the rest of the 
terms in (6.4) and (6.3) arc given solely through the i?-matrix elements /(A) and therefore depend 
only on the related class (2.6) or (2.7). Recall that in integrable models, as described in sect. 3, the 
i2-matrix remains same for all models belonging to a particular class, while the L-operators differ 
and may be obtained through various reductions from the same ancestor Lax operator. 

Therefore taking the Lax operator elements in (6.4) and (6.3) as those from the general Lax 
operator (3.1), one may consider the above eigenvalue and the Bethe equation to be the unifying 
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equations for exact solution of all integrable ultralocal quantum and statistical models constructed 
here. Consequently models like the DNLS, SG, Liouville and the XXZ chain together with the 
6-vertex model, belonging to the trigonometric class (2.6) should share similar eigenvalue relations 
with individual differences appearing only in the form of a(A) and /3(A) coefficients. Thus this deep 
rooted universality in integrable systems helps to solve the eigenvalue problem for the whole class of 
models and for the full hierarchy of their conserved currents in a systematic way. Let us present the 
explicit example of XXZ chain with Lax operator (2.13), defining | > as all spin up state which 
gives a(A) = sin-^(A + a), /3(A) = sin-^ A in Bethe equation (6.3) (with a shift A — > A + ^) resulting 

/ sin(A„ + f) \^_ - sin(A.-Ab + a) 

\sin(A, - f ) j ^J;^ sin(A„ - A;, - a) ' ' 



for a = 1, 2, . . . , m. Similarly (6.4) gives the eigenvalue 



^xxz/^\ ■ N/\ , N -fr sin(Aa - A + f ) . -Ji sin(A — Aa + 3§) 
(A) = sm (A + «)n^^^^^^-3^ + sm Ajg ,^^^_^^^a^ 



(6.6) 



yielding for H^xz = ci, the energy spectrum 



m 1 

= A^(A)'A-i(A) U=o=sina^ _ ^Uin^ ,a. +Ncota. (6.7) 

At the limit a — > 0, sin A — > A, when the i?-matrix along with its associated models reduce to the 
rational class, one can derive the corresponding Bethe ansatz results by taking the rational limit 
of the above equations. For example the relevant equations for the isotropic XXX chain can be 
obtained directly from those for the XXZ chain presented above. Intriguingly the corresponding 
result for the NLS lattice model, which belongs to the same rational class, should also show close 
similarity with that of the XXX chain. 



7 Quantum integrable nonultralocal models 

Though many celebrated classical integrable models like KdV, mKdV, nonlinear cr-model, derivative 
NLS etc. belong to the class of nonultralocal models, successful quantum generalization could be 
made only for handful of them. The reason, as mentioned already, is the violation of the ultralocality 
condition. Recall that this condition helps to transit from local QYBE to its global form and 
consequently establish the integrability for ultralocal systems. Therefore the key equations and the 
related formulation for the integrability theory of the nonultralocal models must be suitably modified. 



7.1 Braided extensions of QYBE 

For understanding first the algebraic structures underlying the nonultralocal systems we have to note 
that the trivial multiplication property (3.4) valid for ultralocal models needs to be generalized here 

as [A (S) B){C (8> D) = ^pBc{A{C B)D) where the braiding iJ'bc takes into account the noncom- 
mutativity of B2,Ci. In spite of such braided extension of the multiplication rule, the associated 
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coproduct structure of the underlying Hopf algebra, crucial for transition to the global QYBE, must 
be preserved. Such a braided extension of the Hopf algebra [49, 50] was implemented in formulating 
the integrability theory of nonultralocal models through an unified approach [51]. The basic idea 
is to complement the commutation rule for the Lax operators at the same site with their braiding 
property at different lattice sites. Note however that in general the braiding may differ widely and 
with arbitrarily varying ranges the picture might become too complicated for explicit description. 
Therefore let us limit first to the nearest-neighbor (NN) type braiding 

L2j+i{lJ.)Z2iLij{X) = Lij{\)L2j+i{n) (7.1) 

assuming that the ultralocality holds starting from the next neighbors. A pictorial description of this 
condition is given in Fig. la). 

The local QYBE at the same time must also be generalized to incorporate the braiding relations, 
such that the transition to its global form becomes possible again. Such braided extension of the 
QYBE (BQYBE) compatible with (7.1) takes the form (see Fig. lb) ) 

Ru{\ - n)Z^^^Lij{X)L2j{fi) = Z^^^L2j{n)Lij{X)Ru{X - m)- (7-2) 

We list below the known nonultralocal integrable models that can be described by the above braided 
equations. Note that the quantum i?-matrix appearing here is the same (2.5) as for the ultralocal sys- 
tems. However the additional braiding matrix Z, unlike the i?-matrix seems to be model-dependent 

and generally independent of the spectral parameter, though similar to the i?-matrix it satisfies the 
YBE like equations and might also become spectral parameter dependent for specific models [51]. 

The next step is the global extension of the BQYBE for the monodromy matrix (2.3) and it is 
not difficult to check that due to the braiding relation (7.1), the form of BQYBE is preserved for 
global matrices hke '^'(A) = U']=i Laj{X) (see Fig. Ic)). However since for the periodic boundary 
condition one imposes LaN+i{X) = Lai (A), the Lax operators Laj{X) for j = 1 and j = N again 
become NN entries and hence modify the equation due to the appearance of an extra Z matrix from 
the braiding relation (7.1), leading finally to the global BQYBE 

Ri2{X - ii)Z^^T^ {X)Z^^T2{ii) = Z^^T2{ii)Z^^T^ (A)i?i2(A - ii). (7.3) 

Though this equation is similar to (2.4), the commutation of the transfer matrices ensuring the 
integrability of the systems through factorization of the trace identity becomes problematic due to 
the presence of Z-matrix. Detail discussion of this problem and the classification of the Z-matrices 
allowing factorization is given in [51]. Investigations of some nonultralocal systems from a different 
angle were done in [52] . It is easy to see that from the corresponding equations for the nonultralocal 
models presented above one can recover the known relations for the ultralocal models by supposing 
the braiding matrix Z = 1 (see also the caption in Fig 1) 
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Fig. 1: Pictorial description of a) braiding relation (7.1), b) local braided QYBE (7.2) for the Lax 
operators Laj{Xa) and c) global braided QYBE for Ta'''\Xa) = Yl'j=i ^aji^a), k < N. Note that putting 
Z = 1, i.e. removing braiding by undoing the crossing of dashed lines in above figures la,b,c) one can recover 
the corresponding pictures for the ultralocal models [1], namely ultralocality condition (1.2), local (2.1) and 
global QYBE (2.4), respectively. 
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7.2 List of quantum integrable nonultralocal models 

Nonultralocal models are mostly nonfundamental systems with infinite dimensional representations 
defined in some Hilbert space. They may correspond to integrable models with spectral parameter 
dependent Lax operator and i2(A)-matrix or may describe only nonultralocal algebras having spectral 
parameterless L-operator and R{X)x^^ao ~^ .R^j^-matrix. Nevertheless the nonultralocal quantum 
models listed below should be described through the same braided relations (7.1,7.2,7.3) or their 
corresponding spectral-less form in a systematic way. Therefore we present only the explicit form of 
their braiding matrix Z and the L operator, indicating the class of i?-matrix they belong to. These 
inputs should be enough to obtain all individual equations and derive the related results. 
I. Systems with spectral parameterless i?-matrix 

1. Current algebra in WZWN model [53] 

The model involves the nonultralocal current algebra 

{Li(x), L2(y)} = |[C, Li(x) - L2{yMx - y) + jC6'{x - y) (7.4) 

with Ci2 = 2Pi2 — 1, where P12 is the permutation operator, L = ^(Jo + <^i) with = d^gg~^, is 
the current and g £ SU (N) the chiral field. Discrctized and quantum version of this algebra may be 
cast as the spectral-free limit of the above braided YBE relations with i?+ as the R-matrix, current 
L as the Lax operator and Z12 = Rq2i as the braiding matrix, which takes the form 

For the details and an interesting quantum group relation of this model the readers are refered to 
the original works [53]. 

2. Coulomb gas picture of CFT [54] 

The Drinfeld-Sokolov linear problem : Qx = C{x)Q describing this system may be given in the 
simplest case by the linear operator C{x) = v(x)a'^ — a'^ with a nonultralocal property due to current- 
like relation {v{x),v(y)} = 5'{x — y). Discretized and quantized forms of the current-like operator 
defined through the commutation relations 

[^fe ''"f ] = ±'''|('^fe,m - Sk+i,i), [vk,vi] = i^{Sk+i,i - 2(5fe,; + 5k,i+i) (7.6) 

■ — 3 ■ + 

construct the corresponding discretized linear operator as = e~^'"k" + Ae^fefi^, , which 

similar to the above case satisfies the spectral-free braided YBE and other relations with i?+ as R 
and Z = g~'^^®<^^ as the braiding matrix. Generalization of this model for SU (N) has also been 
constructed similarly in [54]. 
II. Models with rational i?(A)-matrix 

3. Nonabelian Toda chain [55] 

The Lax operator of the model given by 

LkW = M - ^'^^ ~ ^^=-1 Y Ak = gk9k\ Bk = gk+igk \ 9k e SU{N), (7.7) 
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represents nonultralocal integrable model and solves all braided relations including the BQYBE with 
spectral dependent rational -R(A) = P — ihXI and the braiding matrix Z12 = 1 + ih{e22 ® ei2)vr, 
where P and tt are permutation operators. For further details on this model including its gauge 
relation with an ultralocal model we refer to the original work [55]. 

4. Nonultralocal quantum mapping [56] 

The system is described by the Lax operator L„ = V2„V2n-i, with Vn = A„cr"" + (T+ + ^Vn{i + 
o"^), where the discretized operator Vk = involves nonultralocal algebra (7.6) and yields at the 
continuum limit A ^ the current-like field: iAv{x). This nonultralocal quantum integrable 

model satisfies again integrable braided relations with spcctral-dcpcndcnt rational i?(Ai — A2)-matrix 
similar to the above case but now with a spectral dependent braiding matrix ZyiiX-i) = I — 'x^<^~ ® 
a"*" and Z2i(Ai). For generalization of this model to higher rank groups and other details we refer 
again to the original work [56]. 
III. Models with trigonometric i?(A)-matrix 

5. Quantum mKdV model [57] 

This well known nonultralocal model may be raised to the quantum level with discrete Lax; 
operator 

\ -iA{(»'+)-i Wi^ )■ ' ' 

where Wj^ = e™^ with obeying the nonultralocal relations like (7.6). i?-matrix (2.6) and the 
braiding matrix Z12 = Z21 = 5^2°"'^®'^^ are associated with this nonultralocal integrable system [57]. 
Bethe ansatz solution of quantum mKdV and its generalizations can be found in detail in [58]. It is 
seen easily that one can recover the well known Lax operator of the mKdV field model: U{x,^) = 
^{iv{x)a^ + ^0-2) from (7.8) at the field limit when ^/^Av{x), as Lfe = I + AU{x,$,) + 0{A'^). 

6. Quantum light-cone sine-Gordon model 

It is known that this well known equation : d^_u = 2sin2'u may be represented by the zero 
curvature condition: ^-\J-^- — + [U+, U_] = of the Lax pair U± with U_(x) = ^d-u{x)a^ + 

and similarly for U+(a;). Recently quantum as well as exact lattice versions 
of the nonultralocal Lax operator have been constructed [23], which in particular for U_(x) may be 
given in the form 

4~^'"'^(A) = e^fe-"^"^)'^' + A^ (e-^(f^+'*"^+iV+ + e'(^'^+°"^+iV-) , Vn^ = u^+i - uj. (7.9) 

It may be shown also that (7.9) obeys exactly the above BQYBE and the braiding relation with the 
trigonometric i?-matrix (2.6) and the braiding matrix z[2 ^ = e*"^"^®^"^^ and consequently represent 
a genuine quantum integrable nonultralocal model. 

Some other nonultralocal models known in the literature need introduction of braiding beyond 
NN, basic formulation of which can be found in [50, 51]. Examples of such models having same 
braiding between any two different sites are i.) Integrable model on m,oduli space [59], ii.) Supersym- 
metric models [60, 51], iii. ) Braided algebra [49], iv) NUL extension of YBE [61] etc. Their unified 
description can be found in [51, 62]. 
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7.3 Algebraic Bet he ansatz 

The solution of the eigenvalue problem for intcgrable nonultralocal models by diagonalizing the 
transfer matrix may be formulated through algebraic Bethe ansatz exactly in analogy with the 
ultralocal models, whenever the factorization of the trace problem, as mentioned above, could be 
resolved. The key equation that is to be used for nonultralocal models for finding the commutation 
relations analogous to (6.2) in the ABA scheme should naturally be given by BQYBE (7.3). We 
however skip all details of this ABA formulation for nonultralocal models, which can be found in 
explicit form on the example of the nonultralocal quantum mKdV model in [57, 58] . 

7.4 Open directions in nonultralocal models 

Since some of the nonultralocal models like nonabelian Toda chain, WZWN current algebra, mKdV 
etc. described above can be connected to ultralocal models through operator dependent local gauge 
transformation, it would be challenging to discover similar relation, if any, for the rest of the quantum 
integrable nonultralocal models [23]. 

Other challenging problems undoubtedly are the possible quantum integrable formulation of the 
famous nonultralocal models like nonlinear a-model, complex sine-Gordon model, derivative NLS 
equation etc. through braided YBE. 

As we know, there is a remarkable interconnection between the integrable quantum and statistical 
models. However this connection is discovered until now only for the ultralocal models as we have 
also seen here. Therefore it should be a new direction of study to investigate whether there could 
be any meaningful statistical model corresponding to the integrable nonultralocal models described 
here. 

Another problem worth looking into would be to formulate fundamental nonultralocal models, if 
any, which then could be used possibly for generalizing spin and electron models with nonultralocality. 

Anyway since this vast branch of integrable systems has received significantly insufficient atten- 
tion, we may hope to have many hidden excitements in this area. 

8 Concluding remarks 

Quantum integrable systems can be divided into two broad classes, ultralocal (UL) and nonultralocal 
(NUL). Wc have presented here a brief description of such models with references for further details 
and demonstrated that the models belonging to both these classes can be described systematically 
through a set of algebraic relations signifying integrability of these systems. For UL models these 
relations are the ultralocality condition and the QYBE involving Lax operator L and the R -matrix, 
while for NUL models they are extended to braiding relation and braided QYBE with an additional 
entry of braiding matrix Z. The L operator representing an individual model is naturally model 
dependent and the same seems to be true also for the Z matrix. The i?-matrix on the other hand 
is mainly of two types (elliptic case is not considered here), trigonometric and rational depending 
on the class of models that are associated with q-deformcd and undeformed algebras, respectively. 
This induces a significantly model-independent approach also in the ABA method for solving the 
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eigenvalue problem. For UL systems, the theory of which is more developed, one can go beyond 
and prescribe an unifying algebraic scheme for generating individual Lax operators realized from a 
single ancestor model in a systematic way. It would be a challenge to extend the formulation of 
this scheme also for the NUL models. The integrable statistical vertex models can be related to the 
corresponding quantum models, which as a rule belong to UL systems. Possible systematic extension 
of such relation to NUL systems would be another challenging problem. 
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